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Abstract. 

We consider the existence and stability of the hole, or dark soliton, solution to 
a Ginzburg-Landau perturbation of the defocusing nonlinear Schrodinger equation 
(NLS), and to the nearly real complex Ginzburg-Landau equation (CGL). By using 
dynamical systems techniques, it is shown that the dark soliton can persist as 
either a regular perturbation or a singular perturbation of that which exists for 
the NLS. When considering the stability of the soliton, a major difficulty which 
must be overcome is that eigenvalues may bifurcate out of the continuous spectrum, 
i.e., an edge bifurcation may occur. Since the continuous spectrum for the NLS 
covers the imaginary axis, and since for the CGL it touches the origin, such a 
bifurcation may lead to an unstable wave. An additional important consideration 
is that an edge bifurcation can happen even if there are no eigenvalues embedded in 
**—s | the continuous spectrum. Building on and refining ideas first presented in Kapitula 

Oi , and Sandstede [ p5[ and Kapitula |32[ , we show that when the wave persists as a 

Oi i regular perturbation, at most three eigenvalues will bifurcate out of the continuous 

^^ ■ spectrum. Furthermore, we precisely track these bifurcating eigenvalues, and thus 

Q ' are able to give conditions for which the perturbed wave will be stable. For the NLS 

tjD ' the results are an improvement and refinement of previous work, while the results 

for the CGL are new. The techniques presented are very general and are therefore 
applicable to a much larger class of problems than those considered here. 
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1. Introduction 

The standard model for the propagation of pulses in an ideal defocussing nonlinear 
fiber without loss is the cubic nonlinear Schrodinger equation (NLS) 

^-^xx-^+|^|V = 0, (1.1) 

for igR. It supports the dark soliton solution, which is given by 

$(x) = tanh(V). (1.2) 

If loss is present in the fiber, then the dark soliton will cease to exist. Thus, at a 
minimum amplifiers must be used to compensate for the loss. The effects of linear 
loss in the fiber as well as linear and nonlinear amplification of the wave along the 
fiber will be incorporated into the model. The issues to be discussed in this paper 
are the persistence of the dark soliton under perturbation, and the stability of the 
persisting solution relative to the PDE. In this article, we shall concentrate on these 
issues for a particular perturbation. We emphasize, however, that the methods and 
ideas presented herein are general, and they are applicable to a much larger class of 
problems. Here we will consider a perturbed NLS (PNLS) which is given by 

i<j>t ~ x(f>xx -<P+ H 2 </> = H-d 1( p xx + d 2 <\> + d z \cj)\ 2 (j) + d 4 |0| 4 0), (1-3) 

where e > is small and the other parameters are real and of O(l) in e. The 
nonnegative parameter d\ describes spectral filtering, d 2 describes the linear gain 
(d 2 > 0) or loss (d 2 < 0) due to the fiber, and ^3 and di describe the nonlinear gain 
or loss due to the fiber. 

A related equation is the nearly real complex Ginzburg-Landau equation (CGL) 

<k ~ ^xx -4>+ \<f>\ 2 4> = ie(-di<f> xx + d 2 <p + d 3 \<j>\ 2 (j) + d 4 MV), (1-4) 

where again e > is small and the other parameters are real and of 0(1). The CGL 
governs the nonlinear evolution of perturbations of a simple solution of a basic system 
of partial differential equations at near critical conditions, provided that the basic 



system satisfies some generic conditions (Eckhaus [14 ). The CGL has been proven 
valid in an asymptotic sense for a large class of systems (Collet and Eckman n [|7| , van 
Harten @, Bohemian et al @, Mielke and Schneider Q, and Schneider @f 0). 
The CGL results from an asymptotic expansion, and equation Jl.4| ) with d± = is 
only the 0(1) part of a more extended equation. The inclusion of the d± term is a 
means of modelling the effect of small, nonlinear higher order corrections (Doelman 
0, Popp et al HI, Stiller et al @, §(J) 



It is clear that studying the existence of steady-state solutions to equations (1.3) 
and ( J1.4D amounts to determining the solution structure for the equation 

- I4I' - <f> + |0| 2 = k(i<W + d 2 cb + dM 2 ^ + d 4 \cj>\U) (1.5) 

C = -T-). To do this, one can set 

v ax > ' 



and then study trajectories in the (r, r', ip) phase space. This task has been done in a 




111, Duan et al [13], 
, 53), Marcqet al [38], 



series of papers, of which Doelman and Doelman et al 

Holmes plj , Jones et al pq | , Kapitula and Kapitula et 

and Van Saarloos et al EiJ are a sample. In Section 2 we prove the following theorem 

regarding the persistence of the wave given by (1.2). The result is not entirely new, 

as it is alluded to by Doelman Q . To determine the stability of the perturbed waves 

relative to the PDEs, however, we need more detailed asymptotic information than 



that which is provided in [10 



Theorem 1.1. Suppose that 

da + d 3 + rf 4 = -e 2 cr*(e) - <j, 

where 

2 

9 y " " .") 



a*(0) 



(di + d 3 + ^d 4 ,) 2 (d 1 +d 3 + 2d 4 ). 



Suppose that (e 2 <r*(e) + cr)(di + dj, + 2d±) < 0. If a = 0, then the wave persists as a 
regular perturbation, with the asymptotic expansion 



r{x) 



ip(x) 



${x)+0(e 2 ) 



^((d 1 +d 3 + d 4 )^(x) + ^d 4 ^ 3 (x)) ( (>[<■■'). 



If a 7^ 0, then the wave persists as a singular perturbation. 



Remark 1.2. When a ^ 0, the radial profile of the wave will have a "shelf" 

(SiilH)- 



Remark 1.3. The wave — $, which exists for e = 0, persists under the same conditions; 
our analysis shows that it has the same stability characteristics as $ as well. For 
concreteness, we will simply refer to $ throughout this paper. 

It seems that all previous attempts to consider the stability of the wave, especially 
for the PNLS, have ignored the fact that the wave persists as a singular perturbation 
except on the regular perturbation manifold ^2+^3+^4 = — e 2 cr*; relevant works 
include Burtsev et al jl], Chen et al |J, Ikeda et al [g2[ [2i|], and Lega et al 36 . If the 



parameters do not lie on the regular perturbation manifold, then it may be the case 
that the "shelf" can influence the stability of the wave. One possible way of attacking 
this problem may be through the topological methods first introduced by Jones p4| 
and Alexander et al 111, and later used in a variety of contexts by, for example, Bosc 
et al j|, Doelman et al Q, Gardner and Gardner et al |16[ [17| Q, and Rubin and 
Rubin et al \ 42 , 43] . This issue will not be addressed in this paper and will be a topic 
of future study. 

Here, we suppose that the wave does persist as a regular perturbation. Since the 
equations under consideration are posed on the unbounded real line, the spectrum 
of the linearization about the wave contains continuous spectrum corresponding to 
radiation modes. In addition, the spectrum may contain several isolated eigenvalues 
of finite multiplicity. Because of the translation and rotation invariance of the PNLS 



and CGL, zero is an eigenvalue. It is not, however, an isolated eigenvalue. When 
e = 0, the continuous spectrum for the NLS covers the imaginary axis, while that for 
the CGL covers the negative real axis. Furthermore, there are no point eigenvalues in 
the open right-half plane for either equation. For e ^ 0, the origin is still contained in 
the continuous spectrum. By choosing the parameters appropriately, one can bound 
the continuous spectrum in the closed left-half plane. To determine the stability of the 
wave for e ^ 0, it is thus necessary to locate the point eigenvalues. There are standard 
tools available which can be used to determine the fate of isolated eigenvalues (see, 
for example, Kapitula 132]). However, it is a difficult and nonstandard problem to 
determine the conditions under which eigenvalues can bifurcate out of the continuous 
spectrum, i.e., conditions under which an edge bifurcation can occur. The primary 
issue of this paper is the detection of such eigenvalues. We emphasize that an edge 
bifurcation may occur even if the corresponding eigenfunctions in the unperturbed 
problem are not localized. 

We now turn to an outline of our approach for locating eigenvalues. In many 
respects it follows that presented in Kapitula et al J35|] , which deals with the stability 
of solitary wave solutions for the focusing NLS. The major tool that we use is the 
Evans function, E(X). The Evans function is a complex- valued function depending on 
A £ C with the property that E(X) = whenever A is an isolated eigenvalue. It is only 
defined a-priori away from the continuous spectrum, so it is not immediately clear that 
it can be used to locate embedded eigenvalues and detect edge bifurcations. However, 
as an application of the Gap Lemma, discovered simultaneously and independently 



by Kapitula et al J35[| and Gardner et al 19 1, the Evans function can be analytically 



extended across the continuous spectrum. The analytic extension can then in theory 
be used to locate embedded eigenvalues and to track them under perturbation. 

In the problems considered so far, it turns out that the continuous spectrum 
corresponds to a branch cut for the Evans function. Furthermore, in these problems 
it is only at the branch point that the Evans function has an embedded zero, so only 
from there can an eigenvalue bifurcate. For the problems under consideration both in 
this paper and in Kapitula et al pq | , when e — the edge of the continuous spectrum 
is a branch point of order one, i.e., near the edge of the continuous spectrum we can 
write E(\) = f(\/\ — At,), where /(•) is analytic and \b is the branch point. In |55| 
the stability of the solitary wave to the perturbed focusing NLS was considered. It 
turned out that for a suitably scaled eigenvalue parameter that near the branch point 
At, = iuj the Evans function could be written as 



E(X, e) = \/A -iu + Ae, 

where A £ C depended upon the particular perturbation. Thus, for that problem at 
most one eigenvalue could pop out of the continuous spectrum. 

To determine the location of the zeros of E(X) near Ab for those problems in which 
more than one eigenvalue can pop out of the continuous spectrum, one would like to 
write the Evans function as the scries 



E( 1 ) = J2a nl n , 7 2 = A-A„, 



and then locate its zeros. This task can be accomplished if one can derive asymptotic 
expressions for the coefficients of the series. Fortunately, by suitably modifying the 
ideas and methods of Kapitula |33] , which were developed for doing Taylor expansions 



around isolated eigenvalues, we are able to derive such expressions. Once the zeros 
of the expansion have been located, we take those zeros that lie on the correct sheet 
of the appropriate Riemann surface and invert to find the eigenvalues for the system. 
The interested reader should consult Section 3 for more details. 

It turns out, for both the PNLS and the CGL, that when e — the Evans function 
has a branch point at A = and is nonzero everywhere else in the closed right-half 
plane. Furthermore, when e = the Evans function has the expansion 

E( 1 ) = A 1 a + 0( 1 4 ), 

where 4 g R and 7 is a suitably defined function of A for A near zero (see Section 
3 for details). Thus, for the perturbed problem, there will be three zeros of the 
Evans function near 7 = 0, and hence there will be at most three eigenvalues in this 
region. By computing the lower order terms in the series, we are able to locate these 
eigenvalues and assess the stability of the hole solution. As the following theorem 
illustrates, for the PNLS there are at most two eigenvalues which bifurcate out of the 
branch point A = and leave the continuous spectrum. Furthermore, the d^ term 
must be nonzero (specifically, negative) for the wave to be linearly stable. 

Theorem 1.4- Suppose that d 2 + ^3 + ^4 = — e2(T *( e )j where a* is given in Theorem 
1.1. Also, assume that d 3 + 2g?4 < 0. 



i) Suppose that d\ > 0, and set Pj\ = dj/d\. If 



P31 < —=Pu. - 1, 

5 



then the linearization of fll.3| ) about the perturbed wave yields a positive O(e) 
real eigenvalue given to leading order by 

A- (d \2d\( /l I 4 (l + ^3i+4P 41 /J7 \ 
^--^ + 2d ^[f+g (P 31 + 2P 41 ) 2 'J 

Furthermore, if 

P31 > -f P41 - 1, P31 > -2P<i - 7, 

5 4 

then there is a positive 0(e 3 ) real eigenvalue which is given to leading order by 

7 ,3 



where 



A 2 = ■ e J 

2(P 3 i + 2P 4 i) 



7 = |df(l + P31 + 7 p 4i) 2 (7 + P31 + 2P 4 i) 
9 5 4 



Otherwise, the wave is linearly stable, as no other eigenvalues bifurcate from the 
continuous spectrum (see Figure Q) . 



ii) If d 1 = 0, then the wave is linearly stable as a solution of (1.3) if hd 3 + Ad 4 < 0; 



otherwise, there is an 0(e) eigenvalue which is given to leading order by 
Ai = 




Remark 1.5. The condition that d\ > and d% + 2d 4 < ensures that the continuous 
spectrum is contained in the closed left-half plane for e > and small. 

Remark 1.6. If d 4 = the wave is linearly unstable, with an 0(e) eigenvalue if 
P31 < —1 and an 0(e 3 ) eigenvalue if —1 < P31 < 0. Furthermore, the wave is 
linearly unstable if d 4 > 0. 

Before we discuss the stability of the wave for the CGL, a few comments are in 
order. There have been many recent efforts to determine the stability of the dark 
soliton for the perturbed NLS by using an adiabatic approach (j|, ||, £2| £3|, pq|). 
With the adiabatic approach the wave is predicted to be stable if both d$ + 2d 4 < 
and d\ + d 3 + 6d 4 /5 > hold. If d 4 = 0, then this approach is consistent with the 



result of Theorem 1.4 in that it correctly determines the stability of the wave up to 
0(e). However, it does not predict the existence of the 0(e 3 ) instability; this is not 
surprising, as the adiabatic approach is only meant to understand the dynamics on 
a time scale of 0(l/e). If d 4 ^ 0, then the analysis contradicts the results presented 
in this paper, even at the 0(e) level. This contradiction implies that the original 
adiabatic ansatz for the slow-time variation of the wave is incorrect (see Section 5.5 
for more details). In some way the parameter d 4 has the same effect on the stability 
analysis for the perturbed wave as it has on the solution structure for the steady-state 
problem, i.e., it breaks some kind of "hidden symmetry" (see Doelman pOl). This 
topic would be an interesting avenue for further research. 

When considering the stability of the wave to the CGL, the primary difficulty is 
that the resulting Evans function is not as easy to factor as that associated with the 
PNLS. As such, for general parameter values the location of bifurcating eigenvalues 
cannot be put into an easily readable form. However, one can determine for which 
ranges in the parameter space there will be eigenvalues with positive real part; as with 
the PNLS, it turns out that at most two eigenvalues bifurcate from the continuous 
spectrum. As it can be seen from the following theorem, a primary difference between 
the PNLS and the CGL when considering the stability of the hole solution is the order 
of the eigenvalues. In general, the instability will grow much more slowly for the CGL 
than for the PNLS. 

Theorem 1.7. Suppose that d2 + d% + d 4 = — e 2 cr*(e), where a* is given in Theorem 



1.1. Set 



± _ 3±a-2/3 2 _ 

Mfln ~2 It « ' a ~ 2 

(/*+ = -1.716, fi- n = -1.385). 

i) Suppose that d\ 7^ 0, and set Pj% = dj/d\. If 

d + P 31 + 2P 41 )(1 + P 31 + |p 4 i) < 0, 

2 5 
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then there is one positive real 0(e 4 ) eigenvalue for the linearized problem, and 
the wave is linearly unstable. If 

dl(l + P31 + -P41) > 0, di0t*7n + P 31 + 2P 4i ) > 

5 



di(l + P31 + 3P41) < 0, di(/i+ + P31 + 2P 4 i) < 0, 

5 

then there is a complex pair of 0(e 4 ) eigenvalues with negative real part. 
Otherwise, no eigenvalues bifurcate from the continuous spectrum (see Figure 
g). In either case, if 

(| + P 31 + 2P 41 )(1 + P 31 + |p 41 ) > 0, 
2 5 

then the wave is linearly stable. 

ii) Suppose that d\ = and set 

Q 

0= (d 3 + 2d 4 ) (d 3 + rd 4 ). 
5 

If a > 0, then the zeros of the Evans function inside the curve K are given by 
A 2 ,3 = (-0.595 ± 0.255i) aV, 



and the wave is linearly stable as a solution of (1.4). If a < 0, then the zero of 
the Evans function inside K is given by 

Ai = 1.191 aV, 

and the wave is linearly unstable. 



Remark 1.8. The continuous spectrum remains in the closed left-half plane for all 
values of di, . . . , d& as long as e > is sufficiently small. 



Remark 1.9. The sign of the parameter a corresponds to the manner in which the 
wave is constructed in the (r, r', ip) phase space. The interested reader should consult 
Section 2 for more details. 



Remark 1.10. If di ^ 0, it may be the case that there is a complex pair of eigenvalues 



with negative real part. The interested reader should consult Lemma 4.8 for the 
details. 



The remainder of this paper is organized in the following manner. In Section 2 the 
conditions for the persistence of the wave are derived through the use of dynamical 
systems techniques. In Section 3 we derive the expressions which allow us to compute 
Taylor expansions at the branch point of the Evans function. This section is relatively 
self-contained and can be skipped on a first reading. In Sections 4 and 5 we calculate 
the Taylor expansion for the Evans function for the CGL and the PNLS, respectively. 



Theorem 1.7 follows from Lemmas n. (i and 4.8. Theorem 1.4 follows from Lemma 5.6. 



Section 5 concludes with a brief discussion comparing the approach of this paper with 
the previous adiabatic approaches. 

Remark 1.11. Recently, Li and Promislow |37[ independently and simultaneously used 
some of the ideas present in this paper to study the stability of waves to the equations 
describing pulse propagation in linearly birefringent, lossless fibers. 

2. Existence and persistence 

The steady-state problem for both the PNLS and the CGL is given by 

- \<P" - 4> + MV - it{\di<f>" + d 2 cb + dM 2 ^ + dMU) (2.1) 

(' = ■j-). For existence of the hole solution, which is given by 

$(x) = tanhz (2.2) 

when e = 0, we will want to consider the problem in polar coordinates. Set 

<j>{x)=r{x)e^> {s)ds (2.3) 

to obtain (after dropping higher order terms that do not affect subsequent calculations) 
the three-dimensional system of ODEs 



s' = -2r(l-r 2 ) + rip 2 

-2e 2 d 1 r{d 2 - di + (di + d 3 )r 2 + d 4 r 4 ) 
$ = -2ft/j-2e(d 2 -d 1 + (d 1 +d 3 )r 2 +d A r 4 ). 



(2.4) 



For the system (2.4) there exist are two critical manifolds Aif, which when e = 
are given by 

M$ = {(r, *, V) : r = ±^1 - V 2 A ^ < 2/3}; (2.5) 

we restrict to ip 2 < 2/3 in ([2.5| ) so that the manifolds A4f are normally hyperbolic. 
Each critical manifold of (2.4) has a two-dimensional unstable manifold, W u (A4f), 



and a two-dimensional stable manifold, W s (A4f), which are smooth perturbations of 
the center-stable and center- unstable manifolds which exist when e = |l5|, g5| . As it 
will be seen, it can be shown that W U (M~) fl W s (Mf) ^ 0, and, by the symmetry 
(r,s,ip,x) -> (r,-s,-ip,-x), W U (M+) n W S {M~) ^ 0, both for < e < e for some 
eo > 0. These relationships are clearly satisfied when e = 0, as evidenced by the 
existence of the waves ±$. Assuming that the relevant manifolds intersect, the wave 
$ will persist as long as the parameters are chosen so that critical points exist onM^ 
(also see Doelman f| pf ). Depending how the parameters are chosen, there will be zero, 
two, or four critical points on Mf (counting multiplicities). The condition ip 2 < 2/3 



9 



implies that the critical points onA^^ correspond to stable periodic solutions to fl2.l| ) 

{§ 0- 

To prove the existence of multiple orbits bifurcating from the original heteroclinic 
cycle with the constraint that the orbits remain within an small tube of the original 
cycle, it will be useful to set 

d 2 + d 3 + d 4 = -(e 2 a*+a), (2.6) 

where a* (e) is such that 

<r*(0) - ~l(di + d 3 + U 4 ) 2 (d 4 + d 3 + 2d 4 ), (2.7) 

9 5 



as in the statement of Theorem 1.1. It will henceforth be assumed that the parameter 



er, while small, is independent of e. 

Remark 2.1. Equation ( p.q ) is not a parameter restriction for the CGL, as it can 
always be achieved by going into an appropriate rotating reference frame. However, 
it is a restriction for the PNLS, and determines a balance between the linear loss and 
nonlinear gain. 



Substituting relation (|2j) into the ODE fl2.4|) yields 
r' = s 



s' = -2r(l -r 2 ) + rip 2 

+2e 2 d 1 r[(d 1 + d 3 )(l - r 2 ) + d 4 {l - r 4 ) + e 2 a* + a] 
f = -2^J + 2e[(d 1 + d 3 ){l-r 2 ) + d 4 (l-r 4 ) + e 2 <T*+<r 



2J..IM. _1_ ^„V1 _ ^2\ , j. /-i _ „4\ , c 2„* , J (2-8) 



Since the lowest order at which a appears in (2J5) is at 0(e) in the ^-equation, the 
effect of a on perturbation calculations will only be felt at 0(e + a)e, except in terms 
of the location of critical points on M E , which is discussed below. Hence, for many of 
the perturbation calculations that follow, the role of a can be ignored. 

The following two propositions detail the relevant behavior on M.f. The proofs 
can be found in Kapitula |pl| and hence are omitted. 

Proposition 2.2. Suppose that ^2 + ^3 + ^4 = — (e 2 er* + a) and that 

(e 2 a*+<j){d 1 +d 3 + 2d 4 ) <0. 

Then a pair of critical points on M+ [M~] are given by (r+,0,±ip*) [(r* , 0, ±i/>*)], 
where 



0+ 


1 


e 2 a* 


+ 0- 


2di 


+ d 3 


+ 2d 4 


? 


e 2 o-< 


+ 0- 





V d x + d 3 + 2d 4 



Proposition 2.3. When < e < 1, the manifolds Aif intersect the r-axis. Further, 
there exists 5, with 1 3> 5 > 0, such that for — (ip* + 5) < ip < ip* + S the flow on Aif 
is given by 



1// - e((di +d 3 + 2d 4 )i> 2 + 2e 2 a* + 2a). 
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Proposition 2.2 gives a condition for the existence of critical points on Aif. 
It remains to show that W U (M~) n W S {M + ) ^ for small e ^ 0. Let Eg = 
{(r,s,ip) : r = t/j = 0}. The hole solution belongs to Eg at x = 0, with 
s(0) ^ 0. When e = 0, the manifold W s (Aif) intersects the curve Eg transversely 
in (r, s,-0)-space, since JU s (.Mj) is transverse to the invariant {ip = 0} plane. Thus, 
the intersection will persist for f ^ sufficiently small. Due to invariance under 
(r, s,i/j,x) — > (—r,s,—tp,—x) and the fact that s(0) ^ along the e = solution, it 
can then be concluded that not only does W U (A4~) also intersect Eg transversely, 
but VK U (7W7) n W S (M + ) ^ as well. Hence, the hole solution will persist for e ^ 
and small. The result is not new (for example, see Doelman ||). To determine the 
stability of the wave, however, more information about the wave must be known than 
has previously been given. 

In the remainder of this section, we finish the proof of Theorem (1.1) by showing 
that for a — the perturbed wave arises as a regular perturbation, and then compute 
its asymptotics. We conclude with a discussion of how the nature of the intersection 
that yields the wave differs in various parameter regimes; this is where Proposition 



2.3 is useful. 

Let an underlying hole solution be denoted by (R,S,^). When evaluated at 



e = a = 0, the variational equations associated with (2.8) are given by 



Sr' 


= Ss 


6s' 


= -2(1 - 3i? 2 


H' 


= 2R'V/R 2 5r 




+2[(di 


5e' 


= 


bo' 


= 0. 



-y 2 /2)5r + 2R^5ip 
-2V/Rds-2R'/R8ip 

d 3 )(l-R 2 ) + d 4 (l-R 4 )]Se 



(2.9) 



Since the solution belongs to Eg at x = even for e ^ 0, it is of interest to determine 
the location of the curve Eg as the flow carries it up to the slow manifold A4f. 
Specifically, we wish to determine the -0-coordinates of the points of Eg as they 
approach A4+. Using the fact that the V'-coordinate of Eg is identically zero when 
e = 0, by doing a Taylor expansion we can write that ip — ip c t+0{t 2 ). From evaluation 
of the variational equations over the e = hole solution <&, we find that tp e satisfies 
the initial value problem 

($Ve)' = 2[{d l + ds)(l - 4> 2 ) + d 4 (l - $ 4 )]$ 2 (9 w) 

($Ve)(0) = 0. (2 ' 10) 

Upon integrating, it is seen that 

M%) = I ((di +d 3 + (UMx) + ^4* 3 (z)) • (2-11) 

Let < v < 1 be given, and let T v > be such that 1 — $(T V ) = v. That is, T v 
denotes a time when the curve Eg is within 0(y) of the slow manifold A4 + . Upon 
evaluating the expression for i/) e at T v , it is seen that 

Ve(rv) = J(dl + d8 + |d4) + 0(l/). (2.12) 

A 5 



The following proposition has now been proved. 
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Proposition 2.4- At the time T v such that 1 — <&(T V ) = v, the image of the curve 
Eg under the flow is within an 0{v) distance of the slow manifold M.f, and the 
^-coordinates of points on the image of E^ are given by 

2 8 

V> = i~(di +d 3 + -d 4 ) + 0{v)]t + 0(e + <r)c, 

o 

where 0<6,c<l. 

First suppose that ct = 0. As a consequence of the manner in which a* has 



been chosen, an application of Propositions 2.2 and 2.4 yields that the wave will 



persist as a regular perturbation. This is due t o th e fact that the critical points on 



Adf match the expression given in Proposition 2A. The following lemma gives the 
necessary asymptotics for the perturbed wave. The proof is a standard application of 
perturbation theory, and hence will be left to the interested reader. 

Lemma 2.5. Suppose that a = 0. The perturbed wave then arises as a regular 
perturbation and satisfies 

r = $ + r ee e 2 /2 + 0(e 3 ) 

1> = ^e + 0(e 3 ), 



where 



and 



2 / 3 

ip e (x) = - Ud 1 +d 3 + d 4 )<S>{x) + -d 4 <S> 3 (x) 



r ee (x) = -L[-5(10(di + d 3 ) 2 + 40(di + d 3 )d 4 + 39d 2 4 )<S>(x) 
225 

+8d 4 {5{d 1 + d 3 ) + 8d 4 )$ 3 (x) + 3di$ 5 (x) 
+12d 4 (5(d 1 + d 3 ) + 8<k)x&(x)] 

+^-d 1 [2d 4 ${x) - 3(di + d 3 + 2d4,)x$'(x)]&(x). 

Remark 2.6. Note that 

Lim (2r ee ±^ 2 )(x) = Q. 

x — >±oo 

This fact will be important in later calculations which deal with improper integrals. 
For the rest of this paper, set 

i>+ = lim ij) e (x). (2.13) 

x — >+oo 

Note that by symmetry, lim x _ ) ._ 00 ip e (x) = —i/j£. Upon doing a linear stability analysis 
of the critical points on Aif, one notices the following facts. If 

(di + d 3 + -d 4 )(di + d 3 + 2d 4 ) < 0, (2.14) 

5 
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then the wave will be realized as the intersection of a two-dimensional unstable 
manifold with a two-dimensional stable manifold in the three-dimensional phase space. 
Alternately, if 

{di+d s + -di)(d 1 + d3+2d 4 )>0, (2.15) 

5 

then the wave is realized as the intersection of a one-dimensional unstable manifold 

with a one-dimensional stable manifold in the three-dimensional phase space. In other 



words, if equation (2.14) holds, then the trajectory out of the curve Eg intersects the 
strong stable manifold of the point (r+, 0, e?Pt)'' furthermore, the critical point is an 
attractor on the manifold Mf . This is indicated by Proposition |2.3[ which gives the 
flow o n Mf for j^l <§C 1, and by Proposition |2. 4 If the parameters satisfy equation 



( 2.15 ), then the critical point is a repcllor on the manifold A4^ (see Figure [|). As we 



show in Sections 4 and 5, this structure plays a role when discussing the stability of 
the wave. 

Now suppose that a ^ 0. In this case, the wave arises as a result of a singular 
perturbation, since il) e (T u ) ^ ±-0* at leading order in v. If era* > 0, then the resulting 
wave can be thought of as a concatenation of the solution <f> with solutions tracking 
along close to the slow manifolds M^ . The radial profile of the solution will have a 
"shelf" at the point at which it approaches Aif (see j|, [3, ^2| £3| for a discussion of 
the shelf in the context of the NLS and nonlinear optics). Furthermore, the perturbed 
wave will stay within an Q{e) tube of the original (e = 0) wave <&. Now suppose that 



aa * < . If equation ( 2.14 ) holds, then the wave will stay within an 0(e) tube of $. 
If (2.15) holds, however, then the wave will travel along Aif to a critical point (if it 



exists) outside this tube. 

3. Derivatives at branch points 

Consider the linear operator 

L = Bdl+P(x)d x +N(x), (3.1) 

where B is an invertible n x n matrix whose eigenvalues have nonnegative real part, 
and P(x) and N(x) are smooth n x n matrices satisfying 

lim P(x)=P ± , lim N(x) = N±, 

x — >±oo X — >±oc 

with the approach being exponentially fast. Upon setting Y = [u,u'] T , where 
' = d/dx, the eigenvalue equation Lu — Xu can be rewritten as the first-order system 

Y' = M(A,x)Y, (3.2) 

with 

id 

-B- l {N{x) - Aid) -B~ 1 P(x) 

In this section, we define an Evans function for the operator L. We do this under 
assumptions which imply that at least one of the matrices M±(A) := lim^^ioo M(\, x) 
has a pair of eigenvalues that produce a branch point for the Evans function at a fixed 
value of A. In this context, we develop a technique for differentiating the Evans 
function at this branch point. This method then allows us, in Sections 4 and 5, to 
derive perturbation expansions on a Riemann surface for particular Evans functions 
around branch points. These expansions are crucial in locating eigenvalues for the 
corresponding linear operators. 



M (A, x) 
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3.1. General assumptions and definition of Evans function 



Consider the linear eigenvalue problem (3.2) where M(X, x) g C nx n is smooth in 



x 

for each fixed A and analytic in A for each fixed x. The following assumptions will be 
made on M(X,x). 

Assumption 3.1. The matrix M(X,x) satisfies: 

• \mi x —,± 00 M(X,x) = M±(A), with an exponentially fast approach 

• If Re A > 0, then M±(\) has n eigenvalues with positive real part and n 
eigenvalues with negative real part 

• A pair of eigenvalues for M±(\) are ±yjb(\), where fe(A) is analytic at A = 
with b(0) = and b'(Q) ^ 0, while the other 2n — 2 eigenvalues are analytic at 
A = with nonzero real parts 

• When put into Jordan canonical form, M±(0) has the block 

The second of these assumptions is not necessary, but it holds for the applications 
of interest and we make it to simplify notation. The third and fourth assumptions 
imply that a pair of eigenvalues of M± (A) form a branch point of the Evans function 
at A = 0. Later in this section we will slightly relax the third and fourth ass umptions 
such that this holds for only one of the matrices M±(X) (see Remark |3.6|). Taken 



together, the statements in Assumption 3.1 imply that if M(A, x) is derived from the 



first order system representation of a linear operator L, then {0} € a c (L) and is on 
the edge of the continuous spectrum (see also [£34| |3^]). Finally, we note that while 
it will not be done here, it may be possible to extend the theory to the case where 
M± (0) have several Jordan blocks of the type given above. This could be useful when 
discussing the stability of waves satisfying viscous conservation laws (J19[). 

We now construct the Evans function following the ideas presented in [Q . If A is 
not in the continuous spectrum, then the matrices M±(X) have no eigenvalues with 
zero real part. If each has n eigenvalues with positive real part and n with negative 



real part, then it is possible to define solutions Yj(A, x) to equation (3.2) which are 
analytic in A such that for i = 1, . . . , n 

lim |Yi(A,x)| = 0, Y 1 (A,0)A---AY„(A,0)^0, 

X — ►— OG 

and for i = n + 1, . . . , 2n 

lim |Y<(A,a;)|=0, Y„ +1 (A,0) A • ■ ■ A Y 2n (A,0) + 0. 

x — >+oo 

Following Alexander et al [fil , the Evans function is given by 

S(A)=Y 1 (A,0)A---AY 2 „(A,0). 

If -E(Ao) = 0, then there exists a solution to ( |3.2| ) which decays exponentially fast as 
| a; | — > oo, and hence Ao is an eigenvalue for L. 

If A is in the continuous spectrum, then at least one of the matrices M±(X) 
has an eigenvalue with zero real part, and the above construction breaks down. 
Recently, Kapitula and Sandstede |35| and Gardner and Zumbrun [|19| concurrently 
and independently showed that the Evans function can be analytically extended into 
the essential spectrum via the Gap Lemma. The analyticity of the extension fails 



precisely when Assumption 3.1 holds, as in this case the Evans function has a branch 
point. 
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In many applications, one of which was considered in |3q ], the branch point is 
located on the imaginary axis. Thus, under a perturbation of the wave, it is possible 
for eigenvalues to move out of the branch point and into the right-half of the complex 



plane, leading to an instability. In other words, an edge bifurcation may occur 34 . To 
locate any such bifurcating eigenvalues, our strategy is to do a Taylor expansion for 
the Evans function in the vicinity of the branch point and then to locate the zeros of 
the resulting polynomial; to expand appropriately, we must account for the presence 
of the branch point (|39[). In particular, if a point Ao is a branch point of order k — 1 
for the Evans function, then by setting 7 = (A — Ao) 1 ^ one obtains an expansion 
around the branch point of the form 
00 

One can then find the zeros for £(7) and use the inversion relation A = Ao +j k to find 
the zeros for E(X). The inversion must be done very carefully, however, as the zeros 



of the series (3.3) do not necessarily all correspond to eigenvalues for the linearized 
problem ([3^2;). 

Let K C C be a simple closed curve which encircles the branch point Ao, such that 
no zeros of the Evans function belong to K itself. Furthermore, let K be such that 
it encloses all the possible zeros of E(X) which are contained in the right-half plane. 
The existence of such a curve is guaranteed by a result in Alexander et a l [B , To 
be able to write the Evans function as the infinite series given in equation (BjT), one 
must be able to define the Evans function on a fc-sheeted Riemann surface TZk- The 
surface TZk is constructed in the following manner (|39|], J48|). Let Kq, Ki, . . . , Kf.-i 
be copies of K cut along the nonpositive real axis. Let Sf denote the upper and lower 
edges of the nonpositive real axis regarded as the boundary of Kj , and let 

(A - Ao) 1 /* = |A - A | 1/fc exp[i(arg A - A + 2jn)/k] 

on Kj. Now paste 5$ to S^ , &f to 5% , . .. , S^_ 2 to S^_ l7 and finally S' k _ 1 to 8$ . The 
result is a fc-sheeted Riemann surface TZk, with the sheets coming together at the 
branch point A = Ao. The Gap Lemma (flj|, |35|]) implies that the function E(X) 
extends analytically to the surface TZk , and hence the series is valid. For the zeros of 



the series (^J) to correspond to eigenvalues, they must lie on the correct sheet of the 



Riemann surface. In particular, they must satisfy 

■K 7T 

~k K a] " g7 < I' 

so that they are located on the sheet Kq. Zeros of the series on other sheets correspond 



I < ar S7 < p (3-4) 



to the existence of sol ution s of (3^2) that are not eigenfunctions. 

Under Assumption 3.1, the Evans function will be defined on a 2-sheeted Riemann 
surface. To take into account the fact that a pair of eigenvalues of M±(X) has a branch 
point at A = 0, set 

7 2 = b(X). (3.5) 

By the assumptions on the matrices M±(X), for Re A > there exist solutions 
r fi(\ x), i = 1, ...,n— 1, such that \Yf .( 



Y^.(A, x), i — 1, . . . , n — 1, such that | Y^A, x)\ — > exponentially fast as x — > ±00. 



From the third assumption and equation fl3.5| ), there also exist solutions Y^(7,ar) 
which satisfy 

Urn Yf(j,x)e ± ^=vf( 1 ). (3.6) 

X— >±OQ 
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The vectors vf 1 (7) are analytic in 7 and satisfy 

M ± (7K ± ( 7 )=T7^ ± (7)- (3-7) 



Using the definition of 7 from equation (3.5), the Evans function on the Riemann 
surface is given by 



where 



and 



E(i) = m( 7 , x) (Y s - A Yj A Y+ A Y+)( 7 , x), (3.8) 

Y±( 7 , a = (Y± 1 A---AY± l _ 1 )(7, a: ) 



m(7, x) = exp I — / trM(7, s)ds 

We make a further assumption to allow the possibility of bounded and/or 
exponentially decaying solutions to equation ( |3.2| ) at A = 0; this is not a restriction, 
since we allow k = 0, but simply sets up the notation to handle such solutions. 

Assumption 3.2. The slow solutions satisfy Y~(0,#) = Y+(0, cc). Furthermore, there 
exists a k, with < k < n — 1, such that Y7j(0, cc) = Yt^O, x) for i = 0, . . . , k. 

Remark 3.3. If {0} $■ o- c (L), then fe would be the geometric multiplicity of the 
eigenvalue A = 0. 

The functions ~Y fi ( , y,x) are analytic in A at A = 0; hence, their derivatives with 
respect to 7 are related to derivatives with respect to A by the chain rule, and when 
evaluated at A = 7 = satisfy 



(2m)! 

b'(oy 



ml fl*»Y± (0,*) = i^8TVj 4 (0,x). (3.9) 



The solutions Y^ (7, a;) are not analytic in A at A = 0; however, by the assumptions on 
the eigenvalues of M± (A) they are analytic in 7 (|p9fl). Since Y7(0,a;) = Yj"(0, x), we 



have S(0) = from ( |3.8| ). As a consequence of Assumption |3.2| and equation (3i)), we 



expect that d^ k+1 E(0) ^ with d^E(0) = for < j < 2k. Proving this conjecture 
will be the focus of the next two subsections 

3.2. Derivatives of the slow components 



The definition of the Evans function in (3J5) is based on 2n solutions of equation (3.2 



We can specify a related set of 2n linearly independent solutions {ui, . . . , u 2n } to ( |3.2|) 
at A = 0, which are useful for differentiating components of the Evans function, as 
follows. Set Ui(x) = YT^Ojj;) for i — 1, ...,k. The existence of k independent 
solutions which grow exponentially fast as \x\ — > 00 is guaranteed by a result in 
Gardner and Jones |L7|; let u,-(a;), i = k + 1, . . . ,2k be these solutions. Now set 

u 2k+ i{x) = Yj k+l (0,x), i = l,...,n- k-1 
\in+k-i+i(x) = Yf tk+i (0,x), i = l,...,n-k — l. 
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Finally, set u 2n _i(a;) = Y~(0,:c), and let u 2 „(a;) be chosen so that 

m(0,x)ui(x) A •■■ Au 2 „(x) = 1. (3.10) 

Now, m(Q,x)ui(x) A ■•• A u 2n _i(x) induces a solution u 2 4 n (x) to the adjoint 
equation associated with equation J3.2| ); furthermore, u£ n (x) ■ u 2 „(a;) = 1 (pi, p2[ Q). 
In all of the examples having the branch point structure under consideration of which 
the authors are aware, this particular adjoint solution is bounded above and bounded 
from zero as |a;| — ► oo; hence, this will be an assumption. The theory can be 
appropriately modified if this does not hold true. 

Assumption 3.4- There exist positive constants C\ and C 2 such that the adjoint 
solution uf n (x) satisfies C\ < |u^(a;)| < C 2 for all igR. 

To differentiate the Evans function at 7 = 0, it is necessary to derive an expression 
for <9 7 (Yj" — Y+)(0, x) at some value of x. Set 

Z±( 7 ,x)=Y±( 7 , a; )e ± ^, 
and note that for fixed x, 

d,Yf(0 7 x)=d 1 Zf(0,x). 
Following Kapitula and Sandstede pa], write 

Zf( 1 ,x)=vf( 1 )+Yf(0,x)-vf(0)+w ± ( 1 ,x), (3.11) 

where w (7,2:) is assumed to decay exponentially fast as x — ► ±00 and to satisfy 
w ± (0, x) = 0. This ansatz is valid due to equation (3.6). 

The assumption that &'(0) 7^ implies that we can locally write A = &^ 1 (7 2 ), 
which yields that dX/drf = at 7 = 0. Since M(A, a;) is analytic in A, we then observe 
that (9 7 M(0,cc) = 0. Therefore, it can be readily seen that 

848^(0, x)) = M(0,x)d 1 w ± (0,x) 

(3.12) 
+M(0, x)d 1 vf(0) ± Y±(0, a;). 



The nonhomogeneous term in the above equation decays exponentially fast as 
x — > ±00. This can be seen by noting that as a consequence of equation (3.7), 
9 7 i 
Set 



M ± (0)d^vt(0)=Tv?(0). 



G ± (x)=M{0,x)d 1 vf{0)±Yf(0,x). 
Solving equation ( |3.12 ) with variation of parameters (see p2[) yields 

ri-1 

^±(0,0) - ^ C ±Y^(0,0) + c±Y±(0,0) 



i=\ 



r° 

u*(0) / G ± (x)-uf(x)dx 

J±oo 



2k 



(3.13) 



+u 2 „(0) / G ± (x)-u 2 4 n (x)dx. 

J±oo 



17 

Here uf(x) are solutions to the adjoint equation associated with equation (3.2) 
satisfying uf(x) ■ Uy(x) = &y, and c~ are some constants. As a consequence of the 
manner in which the solutions u^(x) were defined, uf(x) decays exponentially fast as 
|*| — ► oo for i = k+1, . . . , 2k; hence, the improper integrals are valid. The observation 
that 

M(0,x)d J vf(0) ■ uf(x) = -d y vf(0) ■ ^nf(x) 

together with the exp onen tial decay of the adjoint solutions uf(x) simplify the solution 
formula in equation ( p.!3|) to 

Tl-l 

^±(0,0) = ^cfY± i (0,0) + c±Y±(0,0) 



2k 



J2 [^±(0)-uf(0)] 11,(0) 



(3.14) 



i=k+l 



+ [d 7 vf(0) ■ «(±oo) - <(0))] u 2n (0). 

Here we note that since Y^(0, x) — u 2 „_i(a;), YJ-(0, x) ■ uf{x) = for j ^ 2n — 1. 
Therefore, upon an appropriate renaming of the constants one sees that 

n-l 



i=l 



2 k 



+ E [ar(tf-O(0)V(0)] Ui(0) 

i=fe+l 

+ [d 1 (vf-v;)(0)-ut(0)] u 2n (0) 

+ [^7(0) • u£(-oo) - d 7 vf(0) ■ <(+oo)] u 2n (0). 
The following lemma has now almost been proved. 



(3.15) 



Lemma 3.5. Suppose that Assumptions 3J, 3/2, and 3^ hold. The solutions YJ-(7, x) 
then satisfy 

n-l 

9 7 (Y S - - Y+)(0,0) = ^c±Y± i (0,0)+c s Y7(0,0) 






+ [ay«7(0) • <(-oo) - 5 7 «+(0) • <(+oc)] u 2n (0) 



for some constants c { , c s . 



Proof: As a consequence of equation ( |3.11 ), it follows that 

a 7 (Y7 - y+)(o,o) = d>; - «+)(o) + a>- - w +)(o,o), 

where d 1 (w~ — w + )(0,0) is given in equation (3.15). Plugging in the fact that 



2/i 



3>7 - v +)(0) = ]T [^J" - vt){Q) ■ u? (0)] u,(0) 
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therefore yields the result. 



Remark 3.6. If only one of the matrices M±(\), say M_(A), satisfies Assumption 3.1 
i.e., the other matrix, say M + (A), is such that all of its eigenvalues are analytic in A 
at A = 0, then it is only necessary to compute the relevant term <9 7 Yj(0, 0). One can 
then drop the term <9 7 i;+(0) • u^ n (+oo) in the above lemma. 

3.3. Derivatives of the Evans function 

We are now ready to derive expressions for certain derivatives of the Evans function 



with respect to 7 at 7 = 0. Recall Assumption 3.2, which states that there exist k 



solutions at A — to equation (3.2) which decay exponentially as \x\ — > 00. By the 



construction of the system ( p.2j ) it must then be true that for i — 1, . . . , k 

Y%(0,x) = [i> 1 , i ,^ lti \ T , 

where Lipi.i — 0. We assume that although A = is not an isolated eigenvalue of 
finite multiplicity, we can nonetheless find "generalized eigenfunctions" for A = 0. 

Assumption 3. 7. There exist numbers Oi and functions ipj t i, i = 1, . . . , k, j = 1, . . . , Oj, 
such that 

Furthermore, if j > 2, then I^AmC 2 -)! decays exponentially fast as |a;| — > 00. 



Remark 3.8. If A = were an isolated eigenvalue with finite multiplicity, then the 
exponential decay assumption would hold automatically. Otherwise, it is possible for 
the generalized eigenfunctions to either be bounded away from zero or even grow like 
some power of \x\ as \x\ — > 00 (see Section 3.5). 

Set p = J2i=i a ii an d let 

*a s , l (^) = [V'a s , J ,<, l ] T (3.16) 

for i — 1, . . . , k. Following Kapitula J32| it can be shown that d^(Yj t — Yt i )(0, x) = 
for positive integers a < a^, and 

n-l 

d?(Yj ti - Y+ 4 )(0, x) = Yl 4 Y fj(°> x ) + d * Y s(°> *) 
i=i 

(3.17) 

Ik 

+d 2 nU 2n (x) + a t \ Y^ <d\M(0,x)^ ai ^(x),uf(x)>u J (x), 
j=k+i 

for constants cL , d s , and din- In the above, 

+00 



G(x)-H(x)dx. 

-00 



The integrals are valid due to the fact that the adjoint solutions decay exponentially 
fast as \x\ — » 00. 
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Recall the definition of the Evans function given in equation ( |3.8| ). As a 
consequence of the above discussion and equation ( |3.9D , d" l E(0) = for any positive 
integer m < 2p + 1. Upon using relation (3.9), differentiation yields 



d? +l E(0) 



9 7 (Y s --Y+)Aa T (Y7-Y+)A$ 



where 



and 



and 



(2p+l)! 

nlLi(2oi)i 

(2P+ , 1)! ^( Y 7 - Y+) A c) A (Y7 - Y+) A *, 



9 7 (Y 7 - Y+) = 9- (Y^ - Y+ x ) A • • • A ff? (YJ^ - Y+ k ), 
d x (Yj Y+) = 9- (Y 7l - Y+ t ) A • • • A 9-' (Y,- fc - Y+ fc ), 
${x) = m(0, s) (Y 7fc+1 A • • • A Y^ n _ 1 A Y+ A Y+ fc+1 A • • • A Y+„_ 1 



Substituting the result of Lemma 3.5 and equation (3.17) into this expression, one 
obtains the following theorem. 



Theorem 3.9. Suppose that the assumptions leading to Lemma 3.5 hold, and that 
Assumption 3.7 holds. Then derivatives of the Evans function defined from the linear 



operator L satisfy 



d* p+1 E(0) 



(2p+l) 
b'(0)P 



aD 



where 



and 



D 



= d,v;(0) ■ <(-oo) - d 7 vt(0) ■ <(+oo) 

< a A M* oll i, u£ +1 > ■ • • < QxAf* ail i, < fe > 

< d x M* ah>k , u£ +1 > ■■■ < d x M * Qfc , fc , u^ > 



Remark 3.10. A similiar theorem was proved in Kapitula |32| in the case that A = 
is an isolated eigenvalue with finite multiplicity. 



Remark 3.11. Another case that may arise is that &(0) = b'(0) — 0. Since 6(A) is 
analytic, similar expressions for the derivatives of E(pf) at 7 = can be derived via 
the chain rule; the more zero derivatives b(X) has, the more complicated the results. 
Such an example arises in Section 3.5. 
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34- Example: CGL 



Con sider the linearized problem for the CGL (1.4), given in Section 4 in equation 
(4.2). Upon setting e = 0, the matrix Mq(X,x) is given by 



M (A, a) 



2(A 







l + 3$ 2 ) 

2(A 







l + $ 2 ) 



1 

1 







(3.18) 



It is easy to check here that 6(A) = 2A. Following the procedure leading up to equation 
(3.10), choose the solutions to Y' = Mq(0, x)Y to be 



Ul 



[$',0,$",0] T , u 2 



[u\,Q,ul,Qf 



u 3 = [0,$,0,$'] 



'IT 



U 4 



= [o,w 2 ,o, u i] T 



(3.19) 



{u\(x) = x$(x) — 1, u\(x) = $(x) + x$'(x)). The solution u 2 , which grows 
exponentially fast as x — ► ±oo, is chosen so that 






= -1; 



hence, Ui, • • • , 114 satisfies ( 3.10J) . While it is possible to find an explicit expression for 
u 2 , it is not necessary, and hence will not be done. The adjoint solutions satisfying 



5{j are then given by 



Uo 



= [o,t*t.o, 



,2lT 



uf = [0,-$',0,$] 



(3.20) 



Under the normalization Y^(0, x) = U3(a;), a simple calculation reveals that 

»±(7) = [0,±1,0,-7] T (3.21) 

(recall that 7 2 = 2 A in this case). The result of Theorem |3~9| , with ai = 1 and 
$11 = Ui, then implies that 



= a 7 «7(o) • uf (-00) - a 7 «+(o) • uf (+00) = 2, 



and hence 



5^(0) 



12 



16. 



4-oc 



{&y{x)dx 



(3.22) 



The linearized eigenvalue problem when e = can be written as 
L+p = Xp, L-q = Xq, 



where L± are defined in equation (4.4). As such, we can actually say much more 
about the Evans function. First, both operators L± are self-adjoint, so their spectra 
must be real. Furthermore, since L + Q' = and $' has no zeros, an application of 
Sturm-Louiville theory implies that A = is the largest eigenvalue for L + . Similiarly, 
there are no positive eigenvalues for L_. Therefore, the following lemma holds for the 
Evans function. 
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Lemma 3.12. Suppose that e = 0. Set 7- = 2A. For 7 near zero the Evans function 
has the expansion 

£(7) = ^7 3 + 0(7 4 )- 
Furthermore, the Evans function is nonzero for Re 7 > 0. 

Remark 3.13. As a consequence of this lemma, for a perturbed problem it suffices to 
locate the zeros of the Evans function near 7 = to determine the stability of the 
wave. 



3.5. Example: NLS 

Consider the linearized problem for the PNLS ( |l.3| ), given in Section 5 in equation 
(5.1). Upon setting e = 0, the matrix Mq(A, x) is given by 



M Q (X,x) 



10 

1 

2(-l + 3$ 2 ) -2A 

2A 2(-l + $ 2 ) 



(3.23) 



Choose the solutions Y' = Mo(0,x)Y to be those given in equation (3.19), and let 
the adjoint solutions be those given in equation (3.2C). Define 7 by 



7 2 = 2(1 - Vl-A 2 ), 
so that upon taking the principal square root, 



(3.24) 



A 



-7\/4-7 2 - 



Note that 

A = 7 + 0( 7 2 ) 
for 7 sufficiently small, so that 

8_ d_ 

dX d"f 

at (A, 7) = (0,0). Under the normalization Y^(0,x) = 113(2;), a simple calculation 
reveals that 



vf{l) 



■ o [=F7, T\/4 - 7 2 , 7 2 , 7\/4 - 1 2 ] T ■ 



Thus, the result of Lemma 3.5 implies that 



a 7 (Yj - Y+)(0,0) = 2u 4 (0) + ciui(0) + c 3 u 3 (0). 



(3.25) 



(3.26) 



In this exam ple, 6(A) is given in ( 3.24 ), so 6(0) = 0, but b'(0) = as well. As noted 
in R emark [3.10 , this does not in itself r ule o ut use of a modified form of Theorem 
3.9 . Unfortunately, the result of Theorem |3.9| truly cannot be applied here. Since the 
generalized eigenfunctions are given by 



^1,2(2;) = 





$(x) 



^2,2 (x) = - 



x$'(x) + $(x) 
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the assumption that the generalized eigenfunctions decay exponentially fast as \x\ — > 
oo does not hold. Thus, we must construct the desired solutions directly. Using the 
fact that 

{dxYf)' = M d x Yf + d x M a Yf, 
and that Y^(0, x) — Ui(a;), it is not hard to verify that 

5 A Y±(0,a;) = -U4(0)TU3(x). (3.27) 

Thus, upon solving the equation 

(d 2 x Yfy = M Q d 2 x Yf + 2d x M d x Yf 
by variation of parameters, one finds that 

8 2 X (YJ - Y+)(0, 0) = 4u 2 (0) + ci Ul (0). 
Combining this result with equation ( |3.26| ) implies that when e = 0, 

d*E(0) = 3 9 7 (Y s --Y+)Aa2(Y7-Y+)AY+AY+ 



(3.28) 



= -24. 
The following lemma is now almost proved. 



Lemma 3.14- Suppose that e = 0. Set 7 2 = 2(1 — yl — A 2 ). For 7 near zero the Evans 
function has the expansion 

£( 7 ) = -4 7 3 + 0(7 4 )- 
Furthermore, the Evans function is nonzero for Re 7 > except at 7 = 0. 

Proof: It is shown in Chen et al || that the squared Jost solutions of the Zakharov- 
Shabat eigen-equation, i.e., the squared eigenfunctions, form a complete set. In other 
words, bounded eigenfunctions for the linearized problem exist if and only if A € «R 
(or 7 £ iH). Thus, the Evans function is nonzero for Re 7 > 0, and to complete the 
proof we must show that it is nonzero on the set iR\{0}. 

To this end, we will rewrite the eigenvalue problem in such a way as to fully exploit 
the results presented in B . Letting ip = (/>* , the NLS can be rewritten as the system 

i(f>t ~ \<i>xx - 4> + <f> 2 i> = 

-iipt - \i>xx — ip + H j2 = 0. 
Linearizing about the wave <& yields the system 

i<j)t-\<t>xx~^ + '2<5> 2 <t> + <5> 2 ^ = 

-iipt ~ \i>xx ~ i> + * 2 + 2$ 2 V = 0, 
which, upon setting 

(0,^)^(0»e lpt , 
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induces the eigenvalue problem 

\<f/' + (1 - 2<P 2 )cf> - <P 2 tJj = -p4> 

\ip" + (1 - 2$ 2 )?A - ^ 2 cf> = pifi 

{' = d/dx). 

Since 7 G «R if and only if p G R, we will now explicitly construct the Evans 
function for real p. In the usual way, the eigenvalue system 

Y' = M(p, x)Y 

can be constructed. Set 



Z = p+^Ji + p 2 , 

where the principal square root is taken. Note that p£R implies that £ <E R + , and 
that p = implies that £ = 1. The eigenvalues for the asymptotic matrix Mo(£) are 
given by ±£t/(£), ±Ms(£)> where 

and the principal square root is being taken. The corresponding eigenvectors are given 

by 

vf = [1, e, ± M/ , ±£m/] t , «± - [i, -ve, ±/*„ W^F- 

Now, when Re 7 > 0, Imp < 0, so that for Im£ < we need to define the solutions 
Yg and Y - comprising the Evans function so that 

km (Y± A Y±)(£, ^±0^+^/)* = vf A wf. 

This is done so that the definition of the Evans function is consistent with that given 
in equation (|0|). Using the information presented in H, it can readily be checked 
that 

km (Y s - A YJ)(Z,x)e-(»- + »') x = a(Z)b(£) v+ A v+ 

x — *+oo J J 

where 






2 



Thus, we get that 

B(0 = lim (Y s -AY7AY+AY+)(C,x) 

= a(£)6(£) V" A vj A w+ A ^ . 
Since 

v- A W J A vf A u+ = -4?: (1 + e2)2 3 (1 ^ 2) , 

we see that -E(£) ^ for £ £ R + except when £ = 1. As £ = 1 corresponds to p = 0, 
the proof is complete. ■ 
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Remark 3.15. The functions a(£) and &(£) are related to the transmission coefficient 
for the Zakharov-Shabat inverse scattering problem. 



Remark 3.16. As a consequence of Proposition 2.17 in 35 1, the Evans function will 
remain nonzero for e > and I7I sufficiently large. Therefore, for a perturbed problem 
it suffices to locate the zeros of the Evans function near 7 = to determine the stability 
of the wave. 



4. Perturbation calculations at the branch point: CGL 

In the next two sections we will be using the Evans function to locate the eigenvalues 
that bifurcate out of the branch point. To accomplish this task, we will need to 
perform perturbation calculations for the various coefficients of terms in the series 
expansions for the Evans function. Fortunately, the techniques have been developed 
that will enable us to do so. In Kapitula p2[ , a procedure was described which allows 
one to perform these calculations for expansions about an eigenvalue that is isolated 
with finite multiplicity. This assumption is not valid for the systems considered in this 
paper, as we wish to do perturbation calculations around a branch point; however, all 
is not lost. Kapitula and Sandstede [B5[ showed that it is possible to do perturbation 
calculations around a branch point if a transformation is done on the eigenvalue 
parameter so that the branch point does not move under the perturbation. By 
combining and appropriately modifying the approaches of these two works, together 
with the results in Section 3, we are able to do an expansion around the branch point 
in terms of the transformed eigenvalue parameter. Recall the manner in which E("f) 
is defined in equation (3.8). To compute the coefficients in the Taylor expansion for 



E(j), we will need to be able to compute terms such as <9^(Y7 — Yl~)(0,0) for an 
appropriate value of k. The first three subsections are devoted to this task. 
Henceforth, set 



r 



2d 4 , a = Tip2 



(4.1) 



where ?/>+ is specified by ( 2.13 ) and ( 2.11 ). N ote t hat a i s exa ctly the parameter that 
appears on the left hand side of conditions ( 2.14 ) and ( 2.15 ); that is, the sign of a 
is directly related to the structure of the manifolds whose intersection forms the hole 
solution. 



.1. Preliminaries 



After setting <j) — u + iv in equation (IT), let the perturbation of the wave be written 
in the form 



u + iv = (r + (p + iq))e 



f* i>(s) ds 



(this follows the scheme used in Kapitula |27j]). Here r and tp are giv en i n Lemma |2. 5. 
For e ^ 0, the linearized eigenvalue problem derived from equation (1.4), is given, up 
to 0(e 2 ), by 



A 



1 - e 2 d\ ed x 

-edi 1 - e 2 d\ 



1 „ 

= L + eL e + -e L ee , 



(4.2) 



where 



with 



and 



Ln = 



L+ = -di 



L e = 







l-3$ 2 



L--& 



1-$^ 
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(4.3) 
(4.4) 



and 



L, 



-(tpedx - —Ipe) 



-2$ 2 (di 



6$r ee + V>e 




+2di(ip e d x 



1 
-1 



2d 4 $ 2 ) 



(4.5) 




-1 



2$r 



VO 



€ 



i o 

1 



(4.6) 



-4di$ 2 (di + d 3 + 2rf 4 $ 2 



1 




Note that 



L+$' = 0, L_$ = 0. 



In the above, $ is again given by equation (2.2). 

In the standard way, the expansion for the linear operator L given in equations 
(|i~2l)-(fh6l) yields an expansion for the matrix M(X,x), i.e., M = M +M e e+ M £e e 2 /2. 
It is clear that M(X,x) — > M±(\) as x — v ±oo. The branch point for the Evans 
function, A&, is the A value such that the matrices M±(Xb) have an eigenvalue a\. which 
has geometric multiplicity one and algebraic multiplicity two. A routine calculation 
yields the following proposition. 

ll), the branch point of the Evans function is given 



Proposition 

by 


4.1. 


For a given by 




\ 


1 2 4 

= -2 ae - 


Set 


7 = 






= V /2(A-A 6 ). 



For A close to A& the eigenvalues of M± (A) that have geometric multiplicity one and 
algebraic multiplicity two when A = Xj, are given by 



T7' 



■ 4, 



where 



4 



±ae ^ 



26 

When A = A b , the associated eigenvectors are given by 
^ = Tu 4 (0) + ae 2 u 3 (0). 



Remark 4-2. It should be noted that the location of the branch point does not depend 
on which of M±(A) is being discussed. 



.2. Calculations for Y 



Since Yir(A, x) are analytic in an 0(1) neighborhood of the orig in, for fixed x these 
functions have Taylor expansions. Together with Proposition 



4.1 



I, this implies that 

(Y7-Y+)(A b ,0) = (Y7-Y+)(0,0)+d A (Y7-Y+)(0,0)A fc +O(e 8 ).(4.7) 

The behavior of these sol ution s at A = is fairly well understood. As a consequence 
of the derivative formula ( fo.llj ) , 



a A (Y7-Y+)(0,0) 



< dxM{0,x)ux{x),v£{x) > u 2 (0) 
+ CUl (0) + O(e) 



--u 2 (0) + cu 1 (0) + O(e). 



(4.8) 



for some constant c. In addition, since 



Yf{0,x) 



r'{x) 







{rij)){x) 


t4>+ 


r(x) 


r"(x) 





(ri>y(x) 




r'(x) 



(4.9) 



where 



Vm 



lim ip{ x )> 



it is seen that 



(Y7-Y+)(0,0) = 2^ + 




r(0) 


r'(0) 



(4.10) 



Since r(0) = for all e > 0, it is necessarily true that (Y7 — Yt)(0,0) will be a 
multiple of U3(0) for all e > 0, and hence it will not make a contribution in the 
resulting perturbation calculations for the Evans function. Since |Af,| = 0(e 4 ), the 
following lemma has now been proved. 

Lemma 4-3. The difference in the fast solutions satisfies, to leading order, 
d*(Yj - Y+)(A 6 , 0) = 32a 2 u 2 (0) + cuui(0) + c 34 u 3 (0), 

for some constants c^ and C34. Furthermore, 

di(Yj - Y+)(A b , 0) = cyui(0) + c 3j u 3 (0), j = 0, . . . ,3 

for some constants c\j and c 3j -. 
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4-. 3. Calculations for Yf 1 

In this subsection all of the calculations will be performed at 7 = 0, where 

7 2 = 2(A-A 6 ). (4.11) 

As such, the 7 dependence of solutions will be suppressed. Set 

Zf( Xl e)=Yf(x,e)e- a ±*. 
The rescaled variable then satisfies the ODE 

d x Zf(x,e) = (M(x)-a b h id)Zf(x,e), (4.12) 



and the asymptotic matrices are now such that they have the Jordan block 



1 




at 



7 = for all e > 0. Again following the procedure outlined in Kapitula and Sandstede 
PI, set 



Zf(x,e)= V b ± (e)+Yf(x,0)- V b ± (0) + w ± (x,e), 



(4.13) 



where w (x,e) is assumed to decay exponentially fast as x — > ±00 and satisfy 
w (x,0) = 0. Furthermore, w (x, e) should not be a scalar multiple of Ui(cc). The 
vectors 7y±(e) are given in Proposition 4.1. Since 9 e ?7j_(0) = <9 e aj_ = 0, upon recalling 
that M = M + M € e + M ee e 2 /2, it follows that 

d x (d e w ± (x,0)) = M o (x)d e w ± (x,0) + M e (x)Yf{x,0), 



and 



5 B (^«r t (a: > 0)) = M (x)d 2 w ± (x,0) + M (x)d 2 r ] b ± 
+2M e (x)d e w ± {x,0) 
+ (M ££ (x) - d 2 a b ± id)Yf(x,0). 



(4.14) 



(4.15) 



Proposition 4-4- Given the ansatz in equation ( 4.15 ), the relevant solution to ( 4.14 ) 
satisfies 

a e w ± (a;,0) = 0. 

Proof: This follows immediately from the fact that M e (x)Yf(x,0) = 0. ■ 

Upon solving equation ( 4.15J ) with the variation of parameters formulation, and 
using the facts that 






and 



one obtains 



M (x)d^ b ± ■ u? = -d 2 V b ± ■ d x ui 
M €e (x)Yf(x,0) = $(2$r £e +^ 2 )u 3 (0), 



d 2 e {w- - w+)(0, 0) = [d 2 r) b _ ■ uf (-00) - d 2 e rt\ • uf (+00)] u 4 (0) 
<!> 2 {x)(2$(x)r ee (x) + i> 2 (x))dxu 4 {0) + cm(0) 
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for some constant c. A tedious calculation reveals that 

$ 2 {x)(2$(x)r ee (x) + ^ 2 (x))dx = -2d x ^\ 



combined with Proposition 4.1, this yields the following lemma. 
Lemma 4-5. The difference in the slow solutions satisfies 

9 e (Y s --Y+)(0,0) = 0, 
and 

d 2 (Y; - Y+)(0, 0) = -4(~dx + T)^+ u 4 (0) + c 2Ul (0) 
for some constants ci and c 2 . 
Proof: Following the discussion leading up to the lemma, it is seen that 

d 2 (w~ - W +)(0, 0) = -4( Jdi + T)^+ u 4 (0) + cui(0). 



The conclusion now follows from the ansatz given in equation (4.13) and the results 
and 



of Propositions 4.1 



4.4 



4-. 4- Calculations for the Evans function 
Set 



.1 



jdi + r), a = r^ 



r 



where T is specified by ( |4.1[) . In the sequel, all of the evaluations will be performed at 
(7, x, e) = (0, 0, 0), and the constants Ci will be u nkn own (but irrelevant). 



Since cc = d\ , as a consequence of equation ([1; 



d 2 JY 



7^/ 



YJ) = --U2 + C1U1, 



with 



9 7 (Y7~Y+) = 



/ 



Furthermore, as a consequence of Lemma 3.5, 

<9 7 (Y S - - Y+) = 2u 4 + c 2 ui + c 3 u 3 . 

From Lemmas (O and 4.5 one has, respectively, that 
dt{Y- f 

and 



/ 



32a u 2 + C4U1 + C5U3, 



9 e 2 (Y s --Y+) = -4ou 4 



C 6 Ui. 



We are now in position to write down a perturbation expansion for the Evans 
function. In the following, the e-dependence of the Evans function is being implicitly 
assumed. First, 

d * m = ^e 2 (Y s --Y+)Ad £ 4 (Y7-Y+)AY+AY+ 



-6! a a, 



and 



and 



d%E(0) = a 7 (Y s --Y+)Aa e 4 (Y y T-Y+)AY+AYJ. 



= i 4!a2 - 



3 2 d 2 £(0) = d 2 (Y s --Y+)Ad 2 (Y7-Y+)AY+AY+ 
32^ 
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In addition, recall equation ( 3.22 ), which states that 



flPtf(O) = 16. 

Note that all lower derivatives of E are zero. Based on the above expansions, the 
Evans function can be written as 



E( % e) = |( 7 3 + heW - a 2 £ 4 7 + a 2 ~ae% 



(4.16) 



While the zeros of the Evans function can be found analytically, it is difficult to 
analyze the resulting expressions. When di = 0, so that a = a, however, the roots are 
given by 



71 



-1.839 ae 2 , 72,3 = (0.420 ± 0.606i) ae 2 



(4.17) 



Recall that 7 2 = 2(A — A&), where Xb is given in Proposition 4.1. The roots of E(j, e) 
are valid as eigenvalues if and only if Re 7 > 0. This is due to the fact that the sheet 
K Q of IZk corresponds to the principal part of \/2(A — A&). Thus, if a > 0, then 72,3 
represent the valid zeros of the Evans function, while if a < 0, then 71 is the valid 
zero. Upon using the inversion formula A = 7 2 /2 + Ab, one has the following lemma. 

Lemma J^.6. Suppose that d\ = 0. If a > 0, then the zeros of the Evans function inside 
the curve K are given by 

A 2 ,3 = (-0.595 ±0.255i)a 2 e 4 . 

If a < 0, then the zero of the Evans function inside K is given by 

Ai = 1.191 a 2 e 4 . 



Remark 4-7. As a consequence, the linearized operator has an unstable eigenvalue if 
a < 0. 

Now suppose that d\ 7^ 0, and set Pj\ = dj/di. To find the zeros, it is most 
illustrative to do a standard bifurcation analysis. From the definition of a, it follows 
that there is at least one positive real zero if (3/2 + -P31 + 2i-4i)(l + P3i +8P41/5) < 0; 
otherwise, there is at least one negative real zero. In addition, a saddle-node 
bifurcation occurs on the lines 

P31 + 2P 4 i = fi sn , 
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where 

+ 3±a-2/3 2 , .._-. , . 

a4 = ^ 7- ' , " 2 = - 4.18 

2 1 =p a 2 

(^+, = —1.716, /xj n = —1.385). By checking the sign of 7 when d 1 Ei^f,e) = 0, it is 
seen that the zeros created by the saddle-node bifurcation have the opposite sign from 
those described above. 

If ip^ = 0, then a = a = 0, so that the branch point does not move and the zeros 
of the Evans function remain at 7 = 0. For the rest of the discussion, assume that 
ijy£ 7^ 0. If r = 0, then the zeros of the Evans function are given by 7 = and 
7 = ±ae 2 . Upon using the inversion formula A = 7 2 /2 + A&, it is seen that there is an 
eigenvalue at A = 0, and no eigenvalues with positive real part. Thus, it is expected 
that the plane T = will serve as the critical plane for which an edge bifurcation may 
take place. 

Now assume for the rest of the discussion that T ^ 0. Set 

7 = Tip+e 2 y. 
Solving £(7, e) = is then equivalent to solving 

y 3 + y 2 -Hy + Li = 0, fi = ( ~ 

For this equation, a saddle-node bifurcation occurs when /i = a 2 . For < jj, < a 2 , 
there is one real negative zero, and the other two zeros are complex with positive real 
part. For /1 > a 2 , all of the zeros are real, and two are positive while one is negative 
(see Figure ||). 

Using the definition of the variable y and the inversion formula, it is seen that for 
Re7 > 0, 



2 y 



+ \2A 



First suppose that Tipf < 0. To achieve a positive zero for 7, one must then have 
y < 0. Since y 2 + fi > 0, this then implies that there is a real positive eigenvalue A, so 
that the wave is unstable. Now suppose that r^ > 0. One must then look at those 
roots with Rey > 0. If y is real, then it is clear that the resulting eigenvalues A are 
negative. If y = y\ + iy-x is complex with yi > 0, then by checking that 

Re = 2 2 {y 1 +y 2 + n) > 0, 

y yi + yi 

it is seen that the resulting complex pair of eigenvalues has negative real part. The 



picture is summarized in Figure H. Thus, the following lemma holds; Theorem 1.7 



follows from Lemma 4.6 and this result 



Lemma 4-8. Suppose that d\ ^ 0, and set Pj\ = dj/d\. If 

(I + P 31 + 2P 41 )(1 + P 31 + ^P 4 i) < 0, 
2 5 
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then there is one positive real 0(e 4 ) eigenvalue for the linearized problem, and the 
wave is linearly unstable. If 

di(i + P31 + 7P41) > 0, di(/*~ + P31 + 2P 4 i) > 
5 



di(l + P31 + 7P41) < 0, di(A*+ + P31 + 2P 4 i) < 0, 
5 

then there is a complex pair of 0(e 4 ) eigenvalues with negative real part (/ijr Tl are 
defined in equation ( 4.18 )). Otherwise, no eigenvalues bifurcate from the continuous 
spectrum. 



5. Perturbation calculations at the branch point: NLS 

5.1. Preliminaries 

As in the previous section, let the perturbation of the wave be written in the form 

u + iv = {r + (p + iq))e l J« s . 



For e 7^ 0, the linearized eigenvalue problem derived from (1.3) is given up to 0(e 2 ) 

by 



A 



edi -(l-e 2 d 2 ) 

1 - e 2 d\ edi 



L + eL € 



1 



-e 2 L f 



(5.1) 



where the operators Lq, L € , and L ee are specified in equations (|4.3[ )- (tl.6| ). As 
previously, the expansion for the linear operator L given in equations ( |4.2[ )-( [0| ) yields 
an expansion for the matrix M(\, x) with M(\, x) — > M±(\) as x — > ±00. As in ( |4.l| ), 
we set r = d\ + ds + 2d$ and a = T-0+ . 

Proposition 5.1. The branch point of the Evans function is given by 
A 



2(r-di) 



e 3 . 



For A close to A& the eigenvalues of M±(\) which have geometric multiplicity one and 
algebraic multiplicity two when A = A& are given by 



where 



and 



and 



4Ti + A(7)£T7, 



b 1 2 

a ± = ±ae , 



7= v / A 2 -2e(r-d 1 )A + a 2 e 4 , 



A( 7 ) = (T - di)e + V7 2 + (r - di) 2 e 2 - a 2 e 4 . 
When A = Ab, the associated eigenvectors are given by 
7^ = TU4(0) + ae 2 u 3 (0). 
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Remark 5.2. To ensure that A b < 0, it is necessary that 

T-di = d 3 + 2d 4 < 0. 
This condition is consistent with ||, g, [2^, 23 , and it will henceforth be assumed. 

Remark 5.3. Since we are taking the principal square root, note that up to leading 
order A(0) = A b for all e > 0. 

5.2. Calculations for Y f 

As in Section 4.2, we use the Taylor expansions of Y^ (A, x), centered at A = 0, for x 
fixed at the origin. From (p~9|), 

O e Yf(0,x) = (V^TV'e+Wx) +$(x)^(x)u 3 (0), 
so that 

dxM o (0,x)d e Yf(0,x) = 2$(x)(il>e(x) T ^ + )u 2 (0). 
The expression given in equation ( p.27| ) implies that 

M e (x)d x Yf(0,x) = 2^MveWu 2 (0). 

J <P(x) 

Solving the equation 

(S 2 A Y±)' = Af c\ 2 A Y± + M e d x Yf + d x M d e Yf 
by variation of parameters thus gives 

9 2 A (Y7-Y+)(0,0)^2(^y^ _i-i V ^)d* 






/ S(x)*'(a;) cZz ) u 2 (0) + c lUl (0), 

J — OO / 



which upon integrating yields 

5 2 A (Y7 - Y+)(0,0) = |(di + d 3 + ^rf 4 )u 2 (0) + ciui(0). (5.2) 

Evaluating the Taylor expansions for both Y7 — Yf and d\(Yj — Yf ), centered at 
A = 0, and using the fact that A b = 0(e 3 ) from Proposition ^J yield the following 
lemma (to leading order). 

Lemma 5.4- The difference in the fast solutions satisfies 

df(Yj - Y+ )(A b ,0) = 16r&(V e + ) 2 u 2 (0) + ci 4 ui(0) + c 34 u 3 (0), 

where 

4 
b = di + d 3 + -04, 
5 

for some constants C14 and c 34 . Furthermore, 

0f'(Y7 - Y+ )(A b ,0) = c ljUl (0) + c 3 ,u 3 (0), j = 0,...,3 
for some constants cij and c 3 j. In addition 

3 2 A (Y7 - Y+)(A b ,0) = ^6u 2 (0) + Cl ui(0). 



:« 



5.3. Calculations for Yf 1 



The only difference in the results of Proposition pj and Proposition 4.f arises in the 
expression for the branch point A&. Furthermore, since |A&| < 0(e 3 ) in both cases, the 
fact that it changes does not affect the calculations up to 0(e 2 ). Hence, the proof of 
Lemma 4.5 applies here to give the following result. 



Lemma 5.5. The difference in the slow solutions at 7 = satisfies 

9 e (Y s --Y+)(0,0) = 0, 
and 

dl{Y- - Y+)(0, 0) = -4(idi + T)0+ u 4 (0) + csui(O) 
for some constants c\ and C2- 



5.4- Calculations for the Evans function 
Set 

f= 1 -d 1 + r. 

In the sequel, all of the evaluations will be performed at (7, x, e) = (0,0,0), and the 
constants Ci will be unknown (but irrelevant). Recall that <9 7 = d\; using this fact, 
along with equation ( 3.26| ) and Lemmas |5.4| and 5.5, we can differentiate to obtain 
a perturbation expansion for the Evans function. As in the previous section, the 
e-dependence of the Evans function is being implicitly assumed. First, we find 



<9 e 6 £(0) = 



_6! 

2!4! 



^dl{Y- Y+) A dt(Yj Y+) A Y+ A Y+ 



and 



-6!rf6(^+) 3 , 



d!d 7 E(0) = 



_3! 

1!2! 



^^(YJ - Y+) A d 2 n (YJ - Y+) A Y+ A Y+ 



and 



-3!r&v 

3 



2! 



+ 

e ' 



a,&;£(0) = TJIi a,(Y,--Y+)Aa,' T (Y,--Y;)AY,+ AY; 



= -3 2!6 ' 



In addition, recall equation (3.28), which states that 



d*E(0) = -24. 
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All lower derivatives of E are zero, so based on the above expansions, the Evans 
function can be written as 

£(7, e) = -4( 7 3 + hej 2 - j^f &V+e 3 7 - irf6(^+) 3 e 6 ) 

(5.3) 
= -4( 7 + he)( 7 2 - r>+£ 2 7 - ilT 0A+) 3 e 5 ). 

To leading order, the roots are for the Evans function are thus 

7i = —be, 72 = r\A+e 2 , 73 = -\t{^?<?- (5.4) 

These can correspond to true eigenvalues only if Re 7 > 0. First suppose that b < 0, 
so that 71 > 0. From the transformation given in Proposition 5.1, i.e., 



A( 7 ) = (r - di)c + V7 2 + (r-rfi) 2 e 2 -a 2 e 4 , 
we find, to leading order, the positive eigenvalue 



x ^-^- d ^[f + w^- l ) e - (5 - 5) 

Now suppose that ti(j+ > 0, so that 72 > 0, and set 7! — a 2 e 4 = 7e 4 , where 

7 = di(^+) 2 (^ 1 +d 3 + 2d 4 ). 
One obtains, to leading order, the second eigenvalue 

A2 = ~wU-/> (5 - 6) 

which is only positive if 7 > 0. Finally, independent of its sign, 73 is of too high an 
order to correspond to a positive eigenvalue A; hence, it can be ignored. The following 
lemma has now been proved; this also yields Theorem 1.4. 

Lemma 5.6. Let d 3 + 2d^ < 0. Suppose that d\ > 0, and set Pji — djjd\. If 

P31 <—=Pn-l, 

5 



then there is a positive O(e) real eigenvalue given, to leading order, by equation (5.5). 
Furthermore, if 

8 5 

-P31 > —-Pa — 1] -P31 > — 2P41 — -, 
5 4 

then ther e is a positive 0(e ) real eigenvalue which is given, to leading order, by 
equation (5J3). Otherwise, the wave is linearly stable, as no other eigenvalues bifurcate 
from the continuous spectrum (see Figure |lj). If di = 0, then the wave is linearly stable 
i f 5c fe + 4c?4 < 0; otherwise, there is an O(e) eigenvalue which is given by equation 

(H). 
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5.5. Comparison with adiabatic approach 



There have been many recent efforts to determine the stability of the dark soliton for 
the perturbed NLS by using an adiabatic approach ( (J, ||, £2], |2^, |36| ) . Following Lega 
et al Bq] , write the solution to the perturbed NLS as 



= (kR$(k4) +601+6 



2, 



n 



exp[i(qx — fit + qx + 6> )] exp[ze /J" Ve( s ) ds], 
where 

^ = x-ct + a;o, q = kn-C, Cl = --q 2 - (nR) 2 , R 2 = 1 + fc 2 . 

Following the procedure outlined in Appendix C of |36J, and using the requirement 
that d,2 + d-s + d4 = 0(e 2 ) for the dark soliton to persist as a regular perturbation, one 
finds that for the time scale T = et, 

2 6 5 

k T = -it[diC- {di+d 3 )kK- -d 4 kn 3 (l + -k 2 )}(l + k 2 ) 
6 5 d 

KT = [d 3 ( K 2 - 1) + d 4 (n 4 -l) + {d 3 + 2d 4 )k 2 K 2 

+d 4 k k — —diq ]k —nkx- 

2 1 + k z 

A linear stability analysis of the critical point (k,K, c) = (0,1,0) yields the 
eigenvalues 

2 6 

Ai = 2{d 3 + 2d 4 ), A 2 = --(di + d 3 + -d 4 ). 

3 5 

Thus, with this approach the wave is claimed to be stable if both d 3 + 2d 4 < and 
d\ + ^3 + 6^4/5 > hold. If d 4 = 0, then this analysis is consistent with the result 



of Lemma 5.6 in that it correctly determines the stability of the wave up to 0(e). 
However, and this is not surprising, it does not predict the existence of the 0(e 3 ) 
instability. If d 4 + 1 0, then the analysis is consistent with what was found via the 
adiabatic approach in ||, ||, g2[ [23); however, these all contradict the results presented 
in this paper, even at the O(e) level. This contradiction implies that the original 
ansatz for the slow-time variation of the wave in the adiabatic approach is incorrect. 
In some way the parameter d 4 has the same effect on the stability analysis for the 
perturbed wave as it has on the solution structure for the steady-state problem, i.e., 
it breaks some kind of "hidden symmetry" (see Doelman [|10)). As mentioned in the 
Introduction, this would be an interesting topic for future study. 
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Figure 1. Stability regime for NLS {d\ > 0) 
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Figure 2. Stability regime for CGL (di > 0) 
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Figure 3. Projected flow onto {s = 0} (a = (dj + d 3 + 2d 4 )(di + d 3 + 8(4/5) ) 
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Figure 4. Zeros of E(f, e) (di > 0) 



